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We present 7-ray spectra for a set of Type la supernovae models. Our study is based 
on a detailed Monte Carlo transport scheme for both spherical and full 3-D geome- 
tries. Classical and new challenges of the 7 ray astronomy are addressed. We find 
that 7-rays are very suitable to reveal the structure of the envelope and, thus, they 
allow to probe properties of the nuclear burning front and the progenitor, namely 
its central density and global asphericities. The potential problems are discussed 
for the quantitative comparison between theoretical and observed line fluxes during 
the first few months after the explosion. 
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7-ray observations have long been recognized as a potential, valuable tool for 
supernovae research flClayton, Colgate fc Fishman, 1969| ; |Ambwani fc Suther 



Hand, 1988 ; |Chan fc Lingenfelter, 199l|) . Only 7-rays provide a direct link to 



X 

the Ni distribution which hardly depends on details of the physics and on the 
numerical treatment. Different scenarios can be distinguished by line fluxes 
and profiles, the structure of the progenitors can be probed, and the time of 
the explosion can be determined. 7-rays can provide a good determination of 
the 56 Ni production for nearby SNe la because, nowadays, accurate distances 
of nearby galaxies can be obtained by 5Ceph. Moreover, all sky surveys by 
7-rays may provide an unbiased rate of SNe la. As we will discuss below, 
advances in the fields of optical and IR observations and of the theory have 
helped to redefine the goals of SN research by 7-rays. New challenges emerged 
which emphasize their central role as probes for the 3-D structure of SNe la. 

The results presented are based on our gamma-ray codes 1 for spherical 
( |riofiich, Miiller fc Khokhlov , 199~2|) and arbitrary 3-D geometries ( [Hofiich ^ 
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1) A typical Delayed Detonation model (M(Ni)=0.56 Mo, M=Mch) 
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2) HeD model (MW=Q.63 Mo, M=0.8 M(C/O)+0.22 M(fte$ 
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Fig. 1. Comparison of a typical delayed detonation model (DD201) with a he- 
lium triggered Sub-Chandrasekhar model (HeDIO) (from [Hoflich, Wheeler J.C. 



Khokhlov (1998a] )). The dashed lines correspond to the contribution of 5G Ni. 

Lee , 2001). We assume homologous expansion of density and chemical struc- 
tures calculated by spherical and 3-D hydro simulations. All nuclear decay 
lines of 56 Ni and 56 Co are included. Pair production and bound-free opacities 
are taken into account. 



2 Classical and New Challenges 



Classical questions and problems: In parts, this section is based on our 



previous analyses flHoflich, Muller & Khokhlov , 1992; 


H5flich, Wheeler J.C, 


Khokhlov, 1998a|). For further discussions, we also want to refer to 


Burrows et 


al. (1991 


),Kumagai & Nomoto (1997|) and Pinto, Eastman & Rogers (2001). 



Overall, 7 spectra and their evolution is characterized by a turnover from 
a phase dominated by 56 Ni to 56 Co lines (Fig. 1). The time of the explosion 
can be determined by the ratio between the 56 A^i(0.81MeV) and the 56 Co(0.84 
MeV) lines because it varies strongly with time of the explosion but it hardly 
depends on the model. With time, the envelope becomes increasingly trans- 
parent. The spectral evolution depends sensitively on the density and chemical 
structure of the envelope and, thus, provides a valuable tool for the discrim- 
ination of explosion models. For example, sub-Chandrasekhar mass models 
show as a distinguishing feature an outer layer of 56 Ni which reveals itself 
by high 7-ray fluxes and broads already a few days after the explosion (Fig. 
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Fig. 2. Influence of the central density on 7-ray spectra at the example 
of a delayed-detonation model (Model 5p0z22.20 of Hoflich et al. 2001) for 



p c = 2,4&6 x 10 9 g 



cm 



(top to bottom lines). 



1). At late phases (>100 days), the absolute line fluxes and profiles provide a 
direct measure of the total 56 Ni mass and its distribution. For M C h models, 
the high densities close to the center (> 10 9 g cm -3 ) result in the production 
of neutron-rich iron-group isotopes rather than 56 Ni. The size of this central 
region increases with p c . A knowledge of this property is critical to detect be- 
cause systematic variations in p c can produce an offset (up to 0.2 m , pommgucz 



Hoflich, fc Straniero (200~T|) ) in the brightness decline relation flHamuy et al. 



1995| ), a cornerstone of modern cosmology with SNe la (e.g. [Schmidt et al 



( 1998|) ; [Perlmutter et al. (1999Q ). Optical and near IR-spectra do not allow 
to distinguish isotopes. However, the lack of central 56 Ni reveals itself in by 
flat-topped line profiles (Fig. 2). To detect the variation, we need resolutions 
between 20 to 30 which are well within reach for the upcoming INTEGRAL 
mission. 



We want to mention one problem related to the analysis of observed line 
fluxes during the first months after the explosion. An advantage of 7-ray com- 
pared to optical analyses is that the results are insensitive to details of the 
physics or numerical treatment. E.g. the specific opacities do not depend on 
temperature or density. This advantage is somewhat lost along the way when 
comparing the observations with theoretical predictions. Problems are caused 
by the time- and model-dependent line shifts and widths (~ 10, 000km /s), 
the intrinsic response function of the instrument which is non-Gaussian, and 
they are connected to the actual definition used to determine integrated line 
fluxes from synthetic spectra. The actual value may differ by up to a factor of 
2 even if based on a given synthetic spectrum (e.g. [Muller, Hoflich fc Khokhfov 



1991|) ) . Preferable, the same machinery should be employed for both the syn- 



thetic spectra and the observations to derive accurate values or good upper 
limits for the 7 emission. 

New challenges: During the last few years, observational and theoretical 
methods in the field of SN-research advanced significantly well beyond the 
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Fig. 3. Energy deposition by 7-rays at day 1 (left) and 23 (right) based on our full 
3-D MC gamma ray transport ( [Hoflich fc Lee , 200 1|) (upper panels). The diameter 
of the WD is normalized to 100. In the lower panels, the 7 spectra at day 23 are 
given as seen from various G of 90, 60, 30 and 0° (red, blue, green, pink) and (p. 
In the presence of inhomogeneities, typical fluctuations of about 10 % can be seen. 
Their frequency provides a measure of the scale. The explosion model is based on 
the delayed detonation model 5p0z22.20 (2) with p c = 2 x 10 9 g cm~ s assuming 
chemical inhomogeneities according to Khokhlov (2001). 



point imaginable as little as 10 years ago. Back then, SNe la were discovered 
a few days before maximum or later and, in general, light curves were rather 
uncertain and the spectral coverage was poor. Even for nearby SNe la, dis- 
tances to the host galaxy were uncertain by 20 to 30 % . Now, accurate 5-Ceph. 
based distances of nearby galaxies are available (e.g. |5aha et al. ( 1997| ) ) . Com- 
bined with the establishment of the optical brightness decline relation relation 
flHamuy et al., 19951) , SNe la provide a unique tool for cosmology. Robotic, 
systematic SN-search programs (e.g. LOSS, c.f. [Treffers et al. ( 1997|) ) con- 
tinuously increase the sample both at low, intermediate and high red-shifts. 
Often, SNe la are discovered 4 to 5 magnitudes before maximum light ([Rics ; 



et al. , 1999 ; Aldcring, Knop fc Nugent, 2000 ), almost eliminating the selection 
effects by galactic extinction and providing a tight handle on the rise times. 
Detailed observations of optical and infrared spectra and light curves allowed 
sophisticated analyses and test of scenarios, ruling all but out the once popu- 
lar helium triggered detonations, and strongly favoring deflagration or delayed 
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detonation Mch models or, in some instances, the merger scenario (e.g. [Hofiich 



fc Khokhlov (19961 )). Observation of polarization in SNe la has shown that, 
in general, these objects are fairly spherical ( |Wang, et al. , 200 1| ) with the 



noticeable exception of very subluminous SNe la ( [Howell et al., 2001]) . For the 



first time, a direct connection with the progenitors seems to be within reach. 
In particular, there is mounting evidence for a connection between the prop- 
erties of the progenitor, and the physics of the explosion ( [Hofiich, Wheeler, & 



Thielemann, 1998b| ; |Iwamoto et al., 1999| ; Pommguez, Hofiich, fc Straniero' 



2001| ; |Khokhlov, 2001| ). The recent progress may redefine the role of modern 
7-ray astronomy for the field of supernovae, and some of classical goals may 
have been rendered less compelling. 7-rays are particular valuable to measure 
the 3-D structure of SNe la. A comprehensive list of new goals is beyond the 
format of this paper. We want to address two of the areas, namely, 7-rays as 
tools to reveal properties of nuclear burning front and large scale asymmetries. 

Nuclear burning fronts: Within Mqh models, optical and IR LCs and spectra 
can be reproduced by models in which a (slow) deflagration front turns into 
a detonation (e.g. |Khokhlov (1991| )) or, alternatively, a deflagration front is 



rapidly accelerating as in W7 ( [Nomoto, Thielemann, fc Yokoi, 1984] ). However, 
successful models require parameterized descriptions for the propagation of the 
burning front. For a discussion, see Pommguez, Hofiich, fc Straniero (2001|) 
and references therein. The propagation of a detonation front is well under- 
stood but the description of the deflagration and the deflagration to detonation 
transition (DDT) pose problems. On a microscopic scale, a deflagration prop- 
agates due to heat conduction by electrons. Though the laminar flame speed 
in SNe la is well known, the front has been found to be Rayleigh- Taylor (RT) 
unstable, increasing the effective speed of burning ( [Nomoto, Sugimoto fc Me"o] 
1976Q . Recently, significant progress has been made toward a better under- 
standing of the physics of flames. Starting from static WDs, hydrodynamic 
calculations of the deflagration fronts have been performed in 2-D ( |Lisewski| 
et al., 2000|) , and full 3-D ( |Khokhlov, 200 1| ). These calculations demonstrated 



the complicated morphology of the front. |Khokhlov (2001| ) finds that, while 
the expansion of the envelope becomes almost spherical, the inhomogeneous 
chemical structure will fill about 50 to 70% of the star (in mass). The resulting 
chemical inhomogeneities and their scale depends sensitively on the structure 
of the initial WD, i.e. progenitor and the pre-conditioning of the runaway. If 
a DDT occurs at densities needed to reproduce normal-bright SNe la, most 
of the unburned fuel will be consumed during the detonation phase and, by 
enlarge, the chemical inhomogeneities will be eliminated. However, they will 
survive in pure deflagration models. 

Currently, state of the art 3-D calculations are restricted to the regime of linear 
instabilities, i.e. the flamelet regime. The consistent treatment is restricted to 
the early part of the deflagration phase. In these calculations, The initial accel- 
eration of the the deflagration front is followed by a declining rate of burning. 
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Fig. 4. Directional dependence of the spectrum of an ellipsoidal delayed detonation 
model (see Fig. 3) at day 23 and 69 after the explosion. We assumed an axis ratio 
of 0.8 consistent with the optical polarization data for SN1999by. 

This leave a significant fraction of the WD unburned (« 0.5M Q ). The resulting 



structures cannot account for the observations of typical SNe la flKhokhlov 



2001] ). This problem is well known from spherical deflagration models with 
a low deflagration speed (e.g. DF1, [Khokhlov (1991|) ), and it triggered the 
suggestion of a DDT, or, alternatively, showed the need for continuously ac- 
celerating deflagration front ([Nomoto, Thielemann, fc Yokoi, 1984|) . Despite 
the limitations of current models, the chemical clumps will not mix in the 
subsequent phase. To test for possible effects of chemical inhomogeneities pro- 
duced during the deflagration, we have remapped a 3-D chemical structure on 
a spherical explosion model (Fig. 3). The resulting line profiles show fluctua- 
tions of about 10 %. Their frequency provides a measure of the scale of the 
instabilities which depends on the C/O ratio of the progenitor. Weak or absent 
fluctuations in the 7-ray spectra would be a strong indicator for a DDT. 

Large scale asymmetries: One of the open questions is the nature of sublu- 
minous SNe la such as SN 199 lbg ([Filippenko et al., 1992| ; [Leibundgut et 
al., 1993|) . Among them, SN 1999by is one of the best observed SNe la. In 
addition to the studies of optical light curves, detailed polarization spectra of 
the subluminous SN 1999by have been obtained and analyzed ( [Howell et al. 



200 1| ). Whereas 'normal' SNe la tend to show little or no polarization ([Wang 



Wheeler, fc Hoffich, 199"7|) , this supernova was significantly polarized, up to 
0.7%, indicating an overall asphericity of the photosphere of ~ 20%. This re- 
sult suggests that there may be a connection between the observed asphericity 
and the subluminosity in SNe la. Among others, a possible explanations are 
the explosion of a rapidly rotating WD and its effect on the propagation of 
nuclear flames during the explosive phase of burning, or extensive burning of 
carbon just prior to the runaway ( [Hofiich et al. , 2001 ). Polarization in the 
optical measures the asymmetry of the photosphere but, in general, provides 
little information whether it is caused by a global asymmetry in the density 
structure or in the excitation mechanism, i.e. the 56 Ni distribution. To test 
the sensitivity of 7-rays, we have calculated the transport for the subluminous 
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model of ellipsoidal shape (Figs. 4). The flux and line profiles vary as a func- 
tion of inclination (e.g. by ~ 50% at day 23). In combination with optical and 
IR observations, 7-rays can be a key to answer the question on the nature of 
the asymmetry. 
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